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1. Introduction 

In the current area of flavor physics, LHC's 2011 data taking pushes the amount of available 
experimental data in B physics to a new level and expands the window of precision tests of the 
Standard Model. Due to confinement the processes considered at the experiments must involve 
hadronic initial states. The impact on the indirect search of New Physics is therefore limited by the 
uncertainties on long distance effects. A natural way to estimate these non-perturbative hadronic 
contributions within a few percent of theoretical error is given by Lattice QCD. To carefully treat B 
physics on the lattice one has to keep under control simultaneously the finite size effects and, partic- 
ularly, the discretisation effects (the lattice spacing should be smaller than the Compton wavelength 
of the b-quark) induced by the simulation. In practice it is not possible to control both effects in 
one simulation. 

The ALPHA Collaboration has followed a strategy discussed in detail in [1]: in HQET the hard 
degrees of freedom ~ rm, are removed and taken into account by an expansion in the inverse b- 
quark mass m^^. As discussed in those papers and also in earlier work, the benefit is the suppression 
of large discretisation effects which may arise in hadronic quantities when the theory is regularised 
on the lattice. The difficult aspect of that method is that a matching with QCD, which is the field 
theory beUeved to describe the strong interactions, is needed to fix the parameters of the effective 
theory. This step also takes care of removing all UV divergencies appearing in the effective theory. 
In lattice HQET those come as inverse powers of the lattice spacing and thus have to be removed 
non-perturbatively before the continuum limit can be taken. 

The main advantages of this method compared to other existing ones are: i) the theoretical 
soundness of the approach, in particular the existence of a continuum Umit which is in addition 
numerically reachable, ii) the completely non-perturbative treatment at any order in the l/ni], ex- 
pansion, including the matching between QCD and HQET, iii) the self-consistency of the method, 
meaning that HQET is tested rather than assumed as opposed to what other approaches have to do 
for masses around the charm, iv) the numerical cost, which is comparable to that for other setups, 
as the extra computations needed for the small- volume matching between QCD and HQET are very 
cheap compared to large- volume simulations. 

2. Parameters and observables of Heavy Quark Effective Tlieory 

The HQET Lagrangian including terms of order l/m\, reads 

=%QET (x) = ^stat (x) - (Okin ^kin (x) - (O^pin ^spin (x) , (2.1) 
^stat W = Yhix) A) Vh{x) , ^kin(.^) = Whi.x)iy^Yh{x) , ^spin W = Yh{x)(y ''^Yhix) , (2.2) 

with the lowest order (static) Lagrangian ^stat and the first order (l/mb)-corrections ff^in, ^spin, 
giving the kinetic and spin contribution respectively. These are sufficient to compute the b-quark 
mass to subleading order, but in order to compute the pseudo-scalar decay constant we also intro- 
duce the zero momentum projected time component of the heavy-Ught axial vector current, 

A^^^\xo)=Z^^^^a^l^jArix)+c^M\-)], (2.3) 
= WiixhoYsYhix) , aJI) {x) = Yi{x)'^-r5Yi(,Vf- Vf)Mx) , (2.4) 
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Table 1: Presently used large volume ensembles of the Coordinated Lattice Simulations (CLS) consortium. 

where Vf denotes the spatial components (/ = 1,2,3) of the symmetric covariant derivative. The 
effective parameters of HQET that need to be known hence are 

(0 = (mbare , In [Z^^^^] , \ (Okin , ^spin) . (2.5) 

The additional parameter wJbare in the static theory, which does not appear in (2. l)-(2.4), is the 
energy shift which absorbs the l/a divergence of the static energy. If the matching is performed 
in HQET to order 1/nib, it absorbs a 1/a^ term. The non-perturbative determination of these 
parameters has been presented at last years conference [2] and closely follows the general strategy 
originally discussed in [1]. We refer the reader to these papers for any unexplained notation and 
further explanations. What still needs to be emphasised at this point is: 

• light quarks are simulated with two-flavor non-perturbatively improved Wilson fermions, 

• for the static action we use the so-called HYP2 discretisation [3],^ 

• in the course of the non-perturbative (NP) matching to QCD in small volume of extent L\ ~ 
0.5 fm, various heavy quark masses have been simulated relativistically at renormalization 
group invariant (RGI) heavy quark mass M fixed to [4] 

Z^LiMG {4,6,7,9,11,13,15,18,21} , (2.6) 

ranging from slighty above the charm to beyond the bottom quark region. 

To summarize, the parameters of HQET, Q){M,a), are known for the masses given in (2.6) and the 
lattice spacings a corresponding to j3 € {5.2,5.3,5.5} in large volume, c.f. Table 1. 

Now we use the effective parameters to compute the B meson mass wib, the hyperfine mass 
splitting me* — hjb and the pseudo-scalai^ heavy-light meson decay constant /b (/b^ is left for a 
future application). To first order in the 1/mb expansion our main observables are defined by 



mB = mbare + + COkin E^^" + CO^pin £^P'" , (2.7) 

mB.-mB = -^COspin£'P'", (2.8) 
In(a3/VB V^V2) = ln(Z^Q^'^) +ln{a^/^p''^') + bf'am^ 

+ COkinp'^'" + WspinP^P'" + c2 V^*" , (2.9) 



A final analysis with increased statistics will also take into account results for HYPl which are fully compatible. 
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Figure 1: Results for static ground state energies (HYP2) after applying the GEVP method on 2 ensembles. 



while their equivalents at static order are given by choosing CO = (m^'^jlg,In[Z^'^'],ac^^',0,0) com- 



pared to (2.5). The improvement coefficient b^^ is determined to one loop in [5] and the HQET 
energies ^stat^^kin^^spin ^^jj ^j^^ HQET hadronic matrix elements p'^\ p^"" , p'^''^ , p^^" de- 
pend on the set {mji,a), i.e. the simulated pion masses niji and lattice spacings a. They have 
been measured on a subset of ensembles produced within the CLS effort [6] with Nf = 2 flavors of 
0(a)-improved Wilson-Clover fermions, obeying 



m^L > 4.0 , 



270MeV < m„ < 450MeV . 



(2.10) 



For this reason we expect finite volume effects to be negligible. Details of these large volume 
ensembles including the actual size of the statistics used here are summarized in Table 1 . We get 
HQET energies & hadronic matrix elements by solving the generalised eigenvalue problem, GEVP, 

C{t)Vn{t,to)=Xn{t,to)C{to)Vn{t,to), (2.11) 

for anN xN correlator matrix C. The n^'^ state has eigenvalue A„ with eigenvector v„. As has been 
shown in [7], a systematic expansion of C in HQET is given by C{t) = C^''''(f) + Lx<^C'^(0 + 
0(ft)^) in terms of small expansion parameters (Ox °^ 1/nib with X G {kin, spin, A^^)}, depending 
on the actual quantity. The eigenvalue A„ determines the effective energy of the state while the 
eigenvectors enter the expressions for the matrix elements [7]. As a variance reduction technique 
we employ stochastic all-to-all propagators. The heavy-light interpolating quark bilinears used, 

Ok{x) = v7j^(x)7o75 V^/^' (x) , Yh{x): static quark field 

Otix) = Wi^'\x)Yo75Vh{x) , = (1 +a'A/10)^>i(^) (2.12) 

represent different levels of Gaussian smearing [8] for the light quark field Yi^'^K k = I, . . . ,N, 
with APE smeared links [9, 10] in the lattice Laplacian A. Numerical experiments have shown that 
choosing N = 3 with x (a/O.Sfm)^ G {1,4, 10} fixed, gives good results. For further details 
about the general procedure see [11]. Special care has been taken to control the contribution of 
excited states. In our analysis we have chosen a time range to extract the plateaux such that the 
corrections to Ef^^ are small compared to its statistical error; we found this to be the case for 
?>?()> 0.3 fm. Figure 1 shows two ground state energies obtained in this way. An autocoiTclation 
analysis has shown that existing data can be considered decorrelated to a sufficient degree. 

We arrange all data sets for a combined jackknife analysis with 100 estimators and compute 
(2.7)-(2.9) which now depend on {z,mji,a). For each quantity we perform a joint, continuum 
(a — )• 0) and chiral (m^ — )• m^"'' = m^jo = 135 MeV [12]) extrapolation — j^+CL in short form. 
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Figure 2: Le/f.- HQET to 0(l/mb) data and results of ;if+CL extrap. of B meson mass to the physical point 



for all values of z- Right: Dependence m^'^^^ {z) in the continuum and its matching to the physical B meson. 



3. Results 

b-quark mass: We determine wib by combining the HQET parameters and energies using the static 
or the 0(l/mb) expressions in eq. (2.7). In both cases we then use the global fit ansatz 

mB{z,mj^,a) = B{z) + C ■ ml. + C' ■ ml+D ■ a'^ , (3.1) 

where we either set C' = or C = -3f/{l6nf^) computed with /;r = f^o = 130.4 MeV [12] and 
g = 0.51(2) [13]. The data points and resulting values of wib at the physical point (from C' = 0) 
are shown on the left of Figure 2, while its z-dependence is shown on the right. Knowing the latter 
allows to match our computations of the B meson mass to its physical value, 

mB(z,m^''P,0)|^^,^ = mg P , taking nig ^ = 5279.50MeV [12] . (3.2) 

As a result we obtain the dimensionless RGI b-quark mass Zb = L\Mb. Taking the recent estimate of 
Li = 0.405(18) fm [14] and applying the 4-/3-loop running of the coupling/mass in the MS scheme, 
we obtain the b-quark mass to 

'nb{mb)\%%2 =4.21(13)stat(3)«(6)^GeV, [zf = 13.46(21)stat(14).(18),] (3.3) 

'^b(mb)|"2='2^ = 4.23(13)stat(3)«(6)zGeV. [z^^^^ = 13.40(22)stat(16)«(18),] (3.4) 

We explicitly separate the statistical eiTor and the eiTor coming from the scale setting. The last 
uncertainty is mainly due to the quark mass renormalisation constant Zm, entering z in QCD [4], 
i.e. the accuracy to which the values of z could be fixed. Since the difference of the static and HQET 
result of mb is already small, we can conclude that the truncation error of the latter 0(A^/m^) is 
negligible compared to the statistical error. Comparing our new result for example with mb(mb) = 
4.163(16)GeV [15] or 4.19^Q QgGeV [12] shows good agreement given our current uncertainty. 
Besides an increase in our statistics and the inclusion of two new ensembles, we are currently 
improving our estimate of Li which is the dominating source of error in and are reducing the 
error in the lattice spacings. 

With the physical value of Zb our theory is now uniquely defined and we can make further 
predictions to compare with experiment. There are two ways to proceed from here: a) one can 
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Figure 3: J+CL extrapolations to the physical point. Left: Decay constant fg. Right: Spin splitting. 



combine the HQET parameters with matrix elements and/or energies to build new quantities, per- 
form a ;t-|-CL extrapolation and a subsequent interpolation to Zb or b) one once interpolates the 
HQET parameters to get U5{a) = w(zb, a;), use them in new observables and perform just the ;t^-|-CL 
extrapolation. Doing both serves as a consistency check. 

Hyperfine/spin splitting: The leading contribution to this quantity, eq. (2.8), is of 0(1 /mb) and 
vanishes in the static limit. Its intrinsic truncation eiTor thus is 0(A^/nib). To construct this ob- 
servable we use GJspin(a) which has leading lattice artifacts of 0(a). Since this quantity is the one 
most sensitive to systematic errors, we will also profit most from a combined analysis with the ad- 
ditional HYPl data and the improvements mentioned above. Our results at finite a and simulated 
pion masses for the HYP2 action together with the PDG value of 45.78(35) MeV at the physical 
point are shown in the right panel of Fig. 3. 

B meson decay constant: Finally we take eq. (2.9) and rewrite it to extract /b using U5(a), m^^ 
and a. Heavy Meson Chiral Perturbation Theory (HM;t^PT) predicts a dependence 

/b [mn^a] = S [l - I j^ml \x^{m\)\ +C-ml + D-a^, (3.5) 

of /b on m^i, taking g^fj^ as in (3.1). The resulting extrapolation is shown left in Fig. 3 (solid curve). 
If the chiral logarithm is neglected one gets a leading order (LO) fit ansatz with a linear dependence 
on m\ (dashed line). The additional, half-transparent data corresponds to results obtained in the 
pure static theory. In numbers this reads 

HM;^:PT: /B|^;'i2 = 189(6)stat(5)„MeV , LO: /B|^;'i2 = 194(6)stat(5)aMeV, (3.6) 
mAxVT. /B|"3='2' = 172(6)stat(5)„MeV, LO: /bI^'IY = 176(6)stat(5)«MeV. (3.7) 

Including the chiral logarithm of HM;^PT or not, changes the value at the physical point by a small 
amount only. At the moment we take the average of the two extrapolations in HQET to 0(l/mb), 

/Bl"?fJ = 174(11) (2) MeV, (3.8) 

as the central value and include half of the difference as part of the systematic error. Note that 
our estimate of /b is lower than other estimates presented at this conference, see [16] for a general 
review. 
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4. Conclusions 

We have reported on the status of the ALPHA Collaboration's heavy quark project to extract 
relevant B physics quantities from N{ = 2 lattice simulations in the framework of HQET expanded 
to 0(l/mb). The measurement of HQET energies and matrix elements has been done using the 
GEVP approach on ensembles produced by CLS. For the first time we quote results for fn\, and 
/b in the continuum obtained in large volume with three lattice spacings and seven pion masses. 
These values are still preliminary since we are confident to decrease uncertainties further by using 
additional data that has not been taken into account yet. 

Once we have further improved the accuracy of our results, we will also determine /b^, 
hadronic parameters of the S — B mixing or B ^ ;r semileptonic form factors as well as more de- 
tails of the spectrum of hadrons with a b-flavor. As a long term goal we would also like to include 
the dynamical strange and charm sea quark contributions such that the only remaining systematic 
effect comes from the truncation of HQET. 
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